Spin reversal of a rattleback with viscous friction 
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Effective equation of motion of a rattleback is obtained from the basic equation of motion with 
viscous friction depending on slip velocity. This effective equation of motion is used to estimate 
the number of spin reversals and the rattleback shape that causes the maximum number of spin 
reversals. These estimates are compared with numerical simulations based on the basic equation of 
motion. 



I. INTRODUCTION 

A rattleback, also known as a celt or wobblestone, is 
a type of mechanical top with the curious property of 
spin asymmetry. Nowadays, there are different varieties 
of rattlebacks. Fig. [T] shows a Russian rattleback toy 
called stubborn tortoise. When the rattleback is spun 




FIG. 1. A Russian rattleback toy called stubborn tortoise. 



in the clockwise direction which the tortoise has turned 
to, it continues to spin clockwise until it slows to a stop. 
However, when the rattleback is spun in the anticlockwise 
direction, a self induced oscillation occurs, and the spin 
slows down and eventually its direction is reversed. There 
is an inertial asymmetry, because the tortoise's center of 
mass is shifted from the principal axes of the body surface 
ellipsoid. This inertial asymmetry is responsible for spin 
reversal. 

Many analyses and simulations have attempted to ex- 
plain the dynamics of rattlebacks during the last century. 

The first scientific paper on rattlebacks appeared in 
1898 by Walkerfl]. Assuming dissipation-free rolling 
without slipping, he obtained linearized equations of mo- 
tion for contact point coordinate variables and spin vari- 
ables. He analyzed the instability of the spin magnitude 
and direction by studying a characteristic equation and 
showed the relationship between the direction of spin and 
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the oscillation. For a realistic case, the spin n is small; 
therefore, his analysis can explain the behavior of a rat- 
tleback. 

In 1983, Pascal^, using the same assumptions as 
Walker flj, derived effective equations of motion for the 
slowly varying mode by using the method of averaging 
and clarified a rattleback's reversal mechanism. 

During the same period, Markeev[3| obtained results 
similar to that of Pascal Q. He derived two conservation 
laws and provided a comprehensible explanation of a rat- 
tleback's reversal. Their results were extended to second- 
order averaging by Blackowiak, Rand and KaplanQ. 

Moffatt and Tokieda|5( presented a physically trans- 
parent derivation of the effective equations of motion sim- 
ilar to Markeev's derivation. 

In 1986, Bondi@ extended Walker'sQ results to un- 
derstand how spin evolves for a wide range of geometric 
and inertial parameters of the body. 

For a no-slip dissipation-free case, numerical simula- 
tions [7] 8|@ showed that infinite spin reversals occur; 
however, a real rattleback has finite spin reversals be- 
cause of energy loss by slip friction. Thus, it is important 
to analyze the dynamics of rattlebacks with slip friction. 

Furthermore, Karapetian[ToJ discussed the stability of 
rotation of a heavy asymmetric rigid body on a hori- 
zontal plane(celtic stone); however, he did not discuss 
the number of spin reversals. Garcia and HubbardQ 
discussed the limitations of a no-slip case and analyzed 
the effects of dissipation. They derived augmented equa- 
tions of motion incorporateing lumped models for aero- 
dynamic effects, spinning torque and slipping torque due 
to Coulomb friction force of slip velocity. In reality, the 
contact with the horizontal plane is not a point but an 
area, relative angular motion between the surfaces cause 
these spinning torque and slipping torque. These equa- 
tions were solved numerically. Because these equations 
were too complicated, they presented a simplified model 
of spin. The spin model derived by considering energy 
in the spin, oscillation, and dissipation was successful in 
explaining spin dynamics; however, this model was not 
derived from equations of motion. Thus, the effect of the 
slip velocity is not clear. Furthermore, they did not dis- 
cuss the number of spin reversals and its relationship to 
the rattleback shape. 

The observation of actual behavior of rattlebacks leads 
to three questions. 
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First, for the no-slip case, if the initial spin value no 
is small, the rattle oscillation becomes large and spin 
reversal occurs in theory. In contrast, for the slip case, if 
no is small, the rattle oscillation does not increase, and 
spin reversal does not occur. In this critical value 

of the initial spin n c seems to exist, above which the spin 
reversal occurs. What is the value of n c ? 

Second, a real rattleback has a finite number of spin re- 
versals because of energy loss by slip friction. When the 
value of the coefficient of friction is known, how many 
times does the rattleback reverse and how does the num- 
ber of spin reversals depend on the coefficient of friction? 

Third, it seems that the number of spin reversals de- 
pends on the rattleback shape. One rattleback may re- 
verse only once, whereas another may reverse as many 
as three times. Given that the lower 2 surface of rattle- 
back is defined by an ellipsoid fr + |y + = 1, where 
a > b > c = I , ii a b, the oscillations in the x direction 
become large and spin reversal occurs rapidly. After one 
reversal, the oscillations start in the y direction but do 
not become as large, and spin reversal occurs slowly. In 
this situation, the number of spin reversals is small. In 
contrast, if a = b, the rattleback is a disk, and spin rever- 
sals does not occur. Thus, it seems that a critical ratio 
of a and b give the maximum number of spin reversals. 
What is this value? 

In this paper, to answer the above questions, I consider 
the basic equations of motion containing sliding friction 
and obtain linearized equations of motion containing slip 
velocity. I derive the effective equations of motion by ap- 
proximating this velocity. By applying the same averag- 
ing method as that used by PascalQ, differential equa- 
tions of slowly changing variables are obtained. These 
differential equations are used to analyze the relation- 
ship among the number of spin reversals, the coefficient 
of friction, and the rattleback shape. The spin behavior 
obtained from the simulation of the effective equations of 
motion is compared with that obtained from the simula- 
tion of basic equations of motion. 



II. BASIC EQUATIONS OF MOTION FOR A 
RATTLEBACK WITH VISCOUS FRICTION 
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FIG. 2. Body axes x, y, and z with lengths a, b, and c, 
respectively, principal inertia axes x, y, z and angle 8. 



by angle S such that 



cos S — sin 6 \ ( x 
sin 5 cos S J \ y 



In this paper, it is assumed that 6 is small, such as 
O(10 -2 ). The vector at the contact point P, x p , has 
components x, y, and z. u is the upward unit vector at 
point P. 

When the rattleback rotates, point P is near the point 
(0,0, c). When the rattleback is oscillates, the values of 
x, y are |x| < a and \y\ < b, respectively, and z can be 
expanded by a second-order approximation of — < 1 and 



| < 1 as follows: 



z ~ c 1 
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where the parameters p, q, and s are given by 
' cos 2 S sin 2 S x 



p — c 2 



S = c 2 
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& ~ a" 2 
sin 2 S cos 2 S 



q = c 2 cos 5 sin 5 ( ^ I , 
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(1) 
(2) 
(3) 



In this paper, a rattleback is considered to be a uniform 
ellipsoid of mass rn with a smooth lower surface such that 
a > b > c as follows: 



x 2 



b 2 



At the leading order, the following equations are ob- 
tained: 

B P ~(x,y,l), (4) 
m~ (-px-qy, -qx-sy,-l). (5) 

Assuming that 5 is small, the principal moments of iner- 
tia iio, /20 j and /30 are approximated as follows: 



where x,y, and z are the body axes of the rattleback as 
shown in Fig. [2] The distances of the three axes from the 
center of mass are a, b, and c, respectively. The principal 
inertia axes are x, y, and z. The z-axis is downward and 
coincides with 5-axis. The axes x- and y-axes are rotated 



where 



m/i, Z20 — rah, ho - ml 3 , 



b 2 + c 2 
1 l = : 7 1 2 



a 2 + b 2 
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In the following equations, I consider the basic equa- 
tions of motion for the angular momentum L and the 



velocity of the center of mass v g . In section IVI| numer- 
ical simulations of these basic equations of motion are 
performed. The evolution of angular momentum is gov- 
erned by Euler's equation: 



= x p x (Rf + Ru) . 



(6) 



where Ru is the normal reaction at P, and / is the slip 
friction force at P. The dynamics of the center of mass 
is governed by Newton's equation for the center of mass 
velocity v g : 



m-^Vg = (R — mg)u + Rf. 



(7) 



The slip velocity v p , is related to velocities v = x p x ui 



and v g as follows: 



v p = v g - v . 



(8) 



Because /, v p , and j^v p only have components in the 
horizontal direction, using Eqs. ([7]) and ©, the normal 
reaction R is given by 



R = mg 



dt 



v 



(9) 



The Coulomb law is often used to define sliding friction 
as follows: 



\v„\ 



where v p is the slip velocity and fi is the coefficient of 
friction. However, this definition of sliding friction is un- 
defined at v p = 0; therefore, it is difficult to analyze the 
equations of motion. Thus, in this paper, to understand 
the dynamics of a rattleback with friction, I shall ap- 
ply the viscous friction law which states that friction is 
linearly related to v p as follows: 



/ 



-uv p . 



(10) 



The angular momentum L has components in the prin- 
cipal inertia axes as follows 

L = (JxoWi, IwW2, , 

where are the components of the angular velocity u> 
of the rattleback. u> is given by the equation -^u — 
0. In general, because the unit vectors of the principal 
inertia axes have time dependence, the time derivative 
of vector A is given by 

d A 3 A I A\ 

-A=-A + (o,xA), 



where JjA = (jjjA^ei. Then 



at 



nu, 



(11) 



where n = uj ■ u. The dynamical variables of these basic 
equations of motion are the components of the angular 
velocity uii and the slip velocity v p i. 

The equations of motion are obtained by dividing both 
sides of Eqs. ^ and © by 5 and rescaling of the vari- 
ables as follows: 

_ ? L _ u> 
t = «/gt, L — oj = —, 



5 R v p 
If —. v — 
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A = V9V- 



In addition, the length of c is assumed to be unit length, 
i.e., c = 1. 

When a rattleback is turned lightly by hand such that 
the initial value of rotation is ~ 2-7T rad/s, spin reversal 
occurs. Because g ~ 980cm /s 2 , the spin value is n ~ 
0.201. As spin reversal occurs for this spin value or lower 
spin values, I assume that h — ~ O(10~ 2 ). This 
assumption is the same as that for the neighborhood of 
the position of stable equilibrium discussed by Pascal[2|. 
The value of the order offi ~ O(10~ 2 ) is used in the next 
section to approximate the equations of motion. 

In the case of limiting /1 to infinity, these equations 
of motion lead to the no-slip shown by the nu- 

merical simulation performed in section I VII Therefore, 
to consider a case in which spin reversal occurs several 
times, /1 is assumed to be not as small as ft, ~ O(10 2 ) ~ A. 

In the remaining sections, the tilde symbol above the 
variables is omitted. 



III. EFFECTIVE EQUATIONS OF MOTION 

In this section, by applying the same linearization 
method as that of Walkerflj and PascalQ, the effective 
equations of motion for a rattleback are obtained. I set 
c = 1 and n ~ i ~ O(10 -2 ), as discussed in the previous 
section. 

The contact point x p has components (x,y,z), which 
satisfy ^ < 1, M < 1. The second-order terms can 
be neglected. Assuming that S ~ O(10~ 2 ), a ~ O(10), 
and b ~ 0(1), as assumed in the numerical simulations, 
Eqs. ([T]), ((2|), and ((3|) for the parameters p, q, and s are, 
respectively, approximated as 



1 



1 



These equations imply that the second-order terms of 
Vp x , V^Vi V*i x i an d V^y can neglected. In addition, 
the order of the spin value n ~ O(10 -2 ) is used, as dis- 
cussed in the previous section. The second-order terms 
of n, p, s, and q are also neglected. 

For the orders of x and y, when the time derivative 
of x is considered, x is multiplied by the oscillation fac- 
tor. Because a rattleback, unlike a usual top, does not 
oscillate rapidly, x and x are assumed to have the same 
order. In fact, as Eqs. (pp)) and (pm| will show, the 
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oscillations are almost cosfVi^t), and according to Eqs. 

([34]) and ([35]). fcf ~ and v\ ~ respectively. 

When a = 10 and 6 = 3, as assumed in the simulation 
in section IVI1 1/1.2 are i/f ~ 5 and ~ 3; thus, it is 
assumed that ^1,2 ~ O(l). Therefore, the second-order 
terms of 0(x) ~ O(x), and O(y) ~ O(y), such as psxy 
are neglected. 

In this approximation, by using Eqs. (@}, © and (fTTj). 
w = (wi,W2,W3) has the following components, 

u> ~ (qx + sy, ~px - qy, -n), (12) 

and Vq = x p x a; has components 

Do — (pi + qy — ny, + + to, 0). (13) 

Because the second term of the time derivative of vq, 
ll> x Vo, is second order of n, p, q and s, ^i>o are as 
follows: 

d d 

dt v ° * m v ° 

~ (pi + qy — ny, qx + sy + nx, 0), (14) 

where terms such as ny are neglected because h is the 
second-order terms, as will be shown in Eq. (|26j) . 

In Eq. (JSJ) for R, it is seen from Eqs. (J5J and 031 that 
(li v °) ' u IS ^ e secon d-order term; thus, R ~ m. There- 
fore, setting R = m in Eqs. © and 0, the following 
equations are obtained 



Next, I consider an approximation in which v p is ex- 
pressed by x, y, and n. This approximation is crucial for 
simplifying these complicated equations of motion. Eqs. 
flED, (Ha and OB give av„ = -4vn - 4v D . In addi- 



give fiv p = - Tt v - -^p. 
tion, the numerical simulation in section I VII shows that 
lt v ° s"p! therefore, the following approximation is 
obtained: 



1 d 

a dt 



Vq. 



(20) 



Because the components of vo are given by x, y, and n 
from Eq. (|13l) . the dynamical variables reduce to x, y, 
and n as will be shown in Eqs. ([29|) - ([3~T|) . It is easily 
seen from Eqs. (fT4")) and (|2U)) that w p3 is second-order 
term. In addition, the term u: x v p is neglected in the 
time derivative of v p . Thus, 4zV p has components, 



d 9 c ■ ri\ 

Jt Vp ~ 9t Vp ~ (V> w p2>°)- 

Using Eq. (|2ip . 7V p has components 



(21) 



AT P ~ (-Vp2,Vpi,xv P 2 - yvpi). (22) 

Substituting Eqs. ([T7 | -(|2^ ]l into Eq. 01) gives 

Ji(gi + sy) + ra + y + v p2 = 0, (23) 

J 2 (p'x + qy) - ny + x + v pl = 0, (24) 



where 



—L = m(N + Ng + N p ), 

d 
dt 



— v 9 = f- 



(15) 
(16) 



where 

J lS J 1 + l=*±°, J 2 =I 2 + 1 = ^. (25) 
5 5 

For n, the following equation is obtained 



N = X p X —V , 



N g = x p x u, N p = x p x —Vp. 



These equations contain six dynamical variables: x, y, 
n, and the components of v p . Therefore, it is difficult to 
analyze of dynamics of these variables. 

Using the approximation in Eq. (|12[) for u>, and ne- 
glecting the second-order terms of n, p, q and s, the ap- 
proximation of Jj£ has components 

~ (h(qx + sy), -I 2 (px + qy), -hti). (17) 
At the leading order, TVo has components 
Nq\ — -(qx + sy) - nx, 

N 2 ~ px + qy - ny, (18) 
N 03 ^ q(xx - yy) + sxy - p'xy + n(xx + yy), 

where the approximation n(l — p)xx ~ nxx is adopted. 
At the leading order, N g has components 



{-V, x, q(y 2 - x 2 ) + xy(p - s)), 



(19) 



I^n + q(xx — yy) + sxy — p'xy + n(xx + yy) 

+ (q(y 2 - x 2 ) + (p - s)xy) + xv p2 - yv p i = 0. (26) 

Let Vpi,v P 2 be expressed with respect to x, y and n 
from Eq. ([20]). Using Eqs. 03]) and and adopting 
approximations (1 — -j^) ~ 1 and ~ 0(~j]?) ~ from 
Eq. ([20)) , Dpi is given by 



Vpl - ~-(px + qy - ny), 
A* 



- ~{ny + ^-', 
Similarly, v p2 is given by 

v p2 ■ 



-{nx-—). 

/Lt Ji 



(27) 



(28) 



where the approximation qx — (I — s)y ~ — y is adopted. 



By substituting Eqs. (gT]) and (gSJ in Eqs. p3 ]l -([26 ]l , 
and neglecting 0(—)(< 0(qj) terms, such as (Jiq—^)x ~ 
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Jiqx, at the leading order, the following equations are 
obtained: 

Ji(qx + sy) + nx + y + = 0, (29) 
Ji/j, 

X 

hipx + qjj) -ny + x + — — = 0, (30) 
I 3 h + q(xx — yy) + sxy — pxy + n(xx + yy) 

+ {a{y 2 -X 2 ) + {p- S )xy) + ^L-^- = 0. (31) 

flJl fJ,J 2 

From the above expressions, it is observed that the 
main parts are oscillations of x and y, such as JZ+j^y = 
and x+ -j^x = 0. Moffatt and TokiedaQ suggested that 
the terms J\qx, nx, J 2 qy and ny are crucial in creating 
reverse oscillations. The effect of friction is included in 
the terms and 

Jlfl J2fJ> 

Furthermore, to eliminate terms depending on qx and 
qy, the new variables X and Y are defined as follows: 



where 



J : 







TJP 



(32) 



p q 
q s 



In Eq. (|32j) . matrix T is a rotational matrix that diago- 
nalizes the symmetric matrix Q = J~ 1 P~ 1 J~ 1 , and it is 
defined by 



T = 



cos a sm t 
— sin 9 cos ( 



The matrix Q is expressed as 



1 

A 



JL 
Ji 



The rotational angle 9 is given by 
tan 29 



2qVJ^h 



Jis — J 2 p 
2V2 



6ab~O(10- L ) 



(33) 



The eigenvalues of Q, v\ and v 2l correspond to the fre- 
quencies of X and Y, respectively, and they are given as 
follows: 



vl 

For a > b, by using Eqs. H])-© for the parameters p, 
q and s, respectively, eigenvalues v\ and v 2 are given as 
follows: 



a 



a 2 + 6 

' b 2 
b 2 + 6 



6(1 

„2 N 



6(1 



(34) 
(35) 



The equations of motion with respect to X and Y are 
given by 



(D + F + DRF) 



X\ ( 



Y J 



where 



R^TJ^NJT- 1 , 



The matrix i? has components as follows: 

ai a,2 



R 



where m are defined by 



a 3 a 4 



(36) 



a± — mi — n/ci, 02 = -m2 — nfc2> 

03 = — 7713 + nk 3 , 04 = m 4 + nfc 4 . (37) 

In these equations, fcj are given by 

fci = -1^ sin 6> cos 9(f - f- 1 ) ~ a 2 <5 - 0(1), (38) 



fc 2 = v\(j sin 2 6» + Z" 1 cos 2 6>) 
k 3 = uf(fcos 2 9 + f- 1 sm 2 9) 



5b 
5a 
5a 2 



•0(1), (39) 
O(10), (40) 



fc 4 = -v\ sm6 cos 0(/ - /- 1 ) ~ — ~ 0(1), (41) 

and depending on friction parameter /1 are given as 
follows: 



mi = 2^2, m 2 = — — fc 4 , 

^2 ^2M 



m 3 = -— fci, m 4 



J 2 u v\ 



2 fe 



where / is defined by 



f = \hr = 



Ji V b 2 + 6 ' 



(42) 
(43) 

(44) 



Finally, the effective equations of motion are obtained 
as follows: 



X + v{X + aiX + a 2 Y = 0, 
Y + v\Y + a 3 X + a 4 Y = 0, 

hh - {X, Y) (KD 2 -K- SD) 



Here, S and K are 2x2 matrices as follows: 



k 2 — fc 4 



S 



_ / si s 2 

S3 s 4 



(45) 
(46) 

0. (47) 



(48) 
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where Sj are defined by 



5i = v\(ny/ JiJ 2 k 3 ), 



S2 



v\{n\J JiJ 2 ki + — ), 
S3 = v\(ny/ J\J 2 ki - —), 



Si = v\(n\J Ji J 2 k 2 H ). 



From Eqs. (|32"]l and it is found that these rn, 

are approximated as m, ~ ^T". When n ~ O(10~ 2 ) and 
- ~ 0(1O~ 2 ) are considered, a* take a value of the order 

O(10 _1 ); therefore, I will neglect the quantities of 0(a 2 ) 
in the next section. For the no-slip case - (p is sufficiently 
large), 



m, ~ 0. 



(49) 



Thus, it is observed that the effect of friction is included 
as i in the rrii parameters. 



Y(0) = yo, X(0) = and Y(0) = as follows: 



Ci - 


J 1 9 /(,,, 


c 2 - 


- * 2A 2/o, 


dl - 


■ «? 
oa 2 vi 




■ a? 

J ^0 



f 2^iA 



yo, 



8a 2 ^2 



1 . 04 

2 _ *4i^ 



where A = v\ — v 2 . 

To obtain the behavior of n(t), I consider the time 
average of the terms in the rapidly varying mode such as 
e tVit , which are contained in terms, such as X 2 and XX, 
of the equation for h(t). When the variables are given 
by Qx{t) = F x (t)Si(t) and Q 2 (t) = F 2 (t)S 2 {t), (where 
Fi(t) are rapidly varying functions and Sj(t) are slowly 
varying functions such as e - ^*), the time average of Qi 
and Q 2 ((QiQ 2 )) is defined as follows: 



(QiQ 2 ) = Km i / F 1 F S 

T-HX 1 Jq 

Then, neglecting terms 0(a 2 ), 



dt. 



IV. SLOWLY VARYING MODE 

The motion of a rattleback contains a rapid frequency 
rattling mode and a slowly varying amplitude mode. To 
discuss the number of spin reversals, the slowly varying 
mode should be analyzed. In this section, by approxi- 
mately solving the characteristic equation and consider- 
ing the time average of the rapid frequency mode, the 
equations of the slowly varying mode is derived from the 
effective equations of motion which were obtained in the 
previous section. 

To study the mode contained in X(t) and Y(t), I ap- 
proximately solve the characteristic equation given in 
Eqs. (l4"51 and (HFjI) . In this case, n is considered to be 
a constant, because h is a second-order term from Eq. 
(pTLj). Moreover, the quantities of 0(a 2 ) are neglected as 
explained in the previous section. 

Two modes are observed, which correspond to the rat- 
tling motion of the long axial direction e lVlt e~~ t and 
the short axial direction e W2t e~~*. Then, the mode ex- 
pansions of X(t) and Y(t) are given by 



(X 2 



L _Q_- ai t 



(XX) 



-a\e 



(XX) — ( 



1^0 -ait 

2 



(XY) 



a l J2- ai t _ Q 2^| „,2„-q 4 t 



8a 2 



2A 



■Vo e 



<YX), 



(Y 2 ) ~ V -±e~ a *\ (YY) ~ -^a 4 e 
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(YY) ~ - V 3A e ~^\ 

(XY) ~ 0, (XY) ~ 0. 

From these equations, it is observed that (XX) and (YY) 
are 0(a,i). Thus, the main part of h is obtained as follows: 



hh{t) ~ h((XX) - (X 2 )) - k 4 ((YY) - (Y 2 )) 
= - k -±A 2 + k -±B 2 , 



(50) 



where 



.4 



v{ - lx e 2 , B = wVj — lx e 



(51) 



X(t) =ae lXlt +cJe- ,A i t +c 2 e 



-i\*t 



+ c 2 e 



i\% t 



and * indicates a 



where Ai = ^1 + i%r, X 2 = v 2 

complex conjugate. The coefficients c\, c 2 , d\ and d 2 
are obtained from X(t) satisfying the equation of motion 
(Eq. (|4"5"j) ) and initial conditions such as X(0) — xq, 



From the definitions given in Eq. (|51|) . these variables 
satisfy the differential equations by using the approxi- 
mation n r*> as follows: 



A = -^A, B 
2 ' 



a 4 



B. 



(52) 



In the no-slip case (/x is sufficiently large) , from Eq. (|49l) , 
the parameters a\ and a4 are given as a\ = —nk\ and 
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a4 = nki, respectively. Then, substituting these param- 
eters into Eqs. ([50)1 and (I521) gives equations that cor- 
respond to those obtained by Pascal[l|, MarkeevQ and 
Moffatt and TokiedaQ. 

Multiplying Eq. (f50)) by n and using Eqs. (|37| and 
(|5^|) . the following equation is obtained: 

^(N 2 + A 2 + B 2 ) = - mi A 2 ~ m 4 B 2 , (53) 

where N = \flj,n. For no-slip case (rrii — 0), the variable 
Ei = N 2 + A 2 + B 2 corresponds to the sum of energy 
about the rotation TV, and the amplitudes of the long and 
short axial directions ( A and B, respectively) are con- 
served. Moreover, from Eq. (|52)l . it is observed that the 
variable E 2 = A 1 B is conserved, where 7 = In the 
phase spaces of JV, A and B, E\ corresponds to a sphere 
and E2 corresponds to a quasi-hyperbolic cylinder. The 
trajectories of the system are closed curves that inter- 
sect this cylindrical surface and the sphere. Therefore, 
an infinite number of spin reversals is obtained. 

In contrast, for the slip case, it is observed that the 
energy E\ decreases according to the right-hand-side of 
Eq. ([53]) , which relates to friction. By using the approx- 
imation h ~ 0, it is assumed that 7 is almost constant 
and E2 is almost conserved; then, the trajectory of the 
system decreases as the radius of the sphere of E\ de- 
creases. Therefore, a finite number of spin reversals is 
obtained. 



V. NEW PHENOMENON FOR THE SLIP CASE 

For the no-slip case, rattle vibration increases irrespec- 
tive of how small the initial spin no is, and in theory, 
spin reversal occurs. However, for the slip case, when no 
is small, vibration does not increase, and rotation stops. 
Even if rattle vibration occurs, spin reversal does not. 
Thus, a critical value of the spin n c \ may exist under 
which rattle vibration does not increase, and a critical 
value of the spin n C 2 may exist under which rattle vibra- 
tion increases but spin reversal does not occur. There- 
fore, it is assumed that the number of spin reversals is 
finite due to the existence of these critical spin values. In 
this section, I discuss some facts relating to these values. 



A. Critical spin n c i necessary to increase rattle 
vibration 

For spin reversal to occur, rattle vibration must in- 
crease. This vibration increases according to the factor 
e~ ait and e~ Q4 *. For the no-slip case, these factors are 
e +nk lt and e -«fc 4 t From Eqs ([35]) and gj), it is ob- 
served that nk\ and nk^ have the same sign and if one 
mode increases, the other decreases. In contrast, for the 
slip case, mi and (114 are in a\ and 04; hence, a\ > 
and 04 > occur according to the value of n. In this 



case, both modes decrease, and rattle vibration does not 
increase. The spin values n c \± are given by 

_ mi cos 2 6 + f 2 sin 2 9 

cl+ fci J2 n sin 9 cos 9(f — / _1 ) ' 

_ m 4 f 2 cos 2 9 + sin 2 9 

cl ~ k 4 J 2 usin9 cos9(f — ' 

for fci > 0. For the case of k% < 0, these values reverse 
and are given by n c i + = — ^ and n cl _ = Therefore, 
it is observed that rattle vibration decreases for the initial 
spin no in the range n c i_ < no < u c i+. 

B. Critical spin n C 2 necessary to reverse rattleback 

spin 

When the spin is more than n c \, rattle vibration in- 
creases. However, the spin does not necessarily reverse. 
A critical spin value n C 2 may exist over which spin rever- 
sal occurs. 

When starting with no > 0, although it is near the 
value of no, for a while, rattle vibration begins rapidly, 
and spin decreases. Eventually, the rattleback stops spin- 
ning and then reverses direction, and the spin value de- 
creases to n\ < 0. This spin value m is obtained as 
follows. 

Now, consider the case for fci > 0, fc4 > and no > 
0. At first, the modes A and B exponentially increase 
and decrease, respectively: therefore, from Eq. (|5T))) . n 
approximates to 

h n^- k -^A\ (54) 

Integrating this equation from to(n = no > 0) to t\(n — 
n\) gives 

fc, i? r* 1 

h{n x -n ) = -^ A 2 . (55) 

Z Jto 

Furthermore, Eq. (|53j) approximates to 

±-{hn 2 + v 2 A 2 )^- mi v 2 A 2 . (56) 
at 

Integrating this equation gives 

^j\ 2 c-^-(nl-n 2 ), (57) 

where it is considered that A(to) ~ and A(t%) ~ 0. By 
substituting Eq. ([57)1 into Eq. (1551) . n\ is given by 

2mi 

nx = -n + — — . 58) 
fci 

Then it is observed that if n\ < 0, that is, the initial spin 
n > -jjp 1 , one spin reversal occurs, and if n n < 
no spin reversal occurs. Therefore, it is observed that if 
n>d+ = < n o < TT 1 ! ra ttle vibration increases but 
spin reversal does not occur. Similarly, when no < 0, the 
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condition no < — ^p 1 must be satisfied for spin reversal 
to occur. 

Finally, it is observed that there exist critical spins 
n C 2± as follows: 



n C 2+ 



2mi 
~~kT 



2m a 



n c2 - 



When n C 2- < uq < n c i- or n c i + < hq < n c2+ , rattle 
vibration increases but spin reversal does not occur. For 
the case of fci < 0, it is observed that rt c2 _ and n c2+ are 
switched: n c2 _ = and n c2+ = -^ L - 



C. Number of spin reversals n r 

For the no-slip case, because dynamical energy is con- 
served, an infinite number of spin reversals occur. In con- 
trast, for the slip case, because friction exist, the number 
of spin reversals n r is finite. 

As discussed in the previous section, for n\ < 0, the 
condition for the first spin reversal occurrence is given by 



1 

hi 



where 



hi 



2mi 



kin 



(59) 



(60) 



Next, for the second spin reversal, if rii satisfy rii > 
?i cl _ = — -2—, rattle vibration does not increase and the 
rattleback stops spinning. If ni satisfy n c i_ > ni > 
n C 2- = — , rattle vibration increases but spin reversal 
does not occur. Therefore, the condition for the second 
spin reversal to occur is as follows: 



ni < — - 



(61) 



the condition that spin reversal occurs only two times 
is given by 1 + p < < 2 + p, only three is given by 
2 + p < j~ < 2 + 2p, and so on. 

It is observed that the number of spin reversals n r 
satisfies the following inequalities: 

iljn > 0, ki > 0, k t > or n < 0, k x < 0, fc 4 < 

K (n r -l,p)<^< K a (n r ,p), (63) 
rii 



illjn > 0, fci < 0, fc 4 < or no < 0, fci > 0, fc 4 > 

Ki(n r - l,p) < ^- < Ki(n r ,p), 
rii 



(64) 



where 



K (n,p) = J2P+(i) + P-(i)P^ 

n 

Ki{n,p) = Y J P+{i)p + P-{i)> 



P±(i) 



i±(-iy 



Finally, from the inequalities stated in Eqs. (|63l) and 
(|64p. the number of spin reversals n r are obtained. 
For case (I), n r is given by 



l 

i+p 



2[x] + 1 if x - 1 < [x] < x 
2[x] ii x- < [x] < x 

and for case (II), n r is given by 

2[x] + l ifx-l<[a;]<a;-l + I ^r 
2[x] if x - 1 + < [x] < x 

where x is defined by 



(65) 



(66) 



Now, I shall estimate spin n 2 after the second rattle 
vibration. Similar to the derivation of n 2 in the previous 
section, the following equation is obtained: 



n 2 



-m 



= n 



2m i 
ki 
2mi 



2m 4 
kd 



Thus, if n 2 > 0, the second spin reversal occurs. This 
condition is consistent with Eq. (|6ip . The condition for 
the second spin reversal occurrence is given by 



1 



>1 + P, 



(62) 



where 



hi 
hi 



2to 4 



k^nQ 



From Eqs. (15^1) and (p2^) . the condition that spin reversal 
occurs only once is given by 1 < ^- < 1 + p. Similarly, 



and [x] is Gauss' symbol, which represents the greatest 
integer less than or equal to x. 



D. Relationship between number of spin reversals 
and rattleback shape 

In this subsection, the relationship among the number 
of spin reversals, the axis inclination (5 ~ 9 in Eq. (|33[) h 
the coefficient of friction p and the rattleback shape are 
obtained. 



The factor x = 



Mp+i) 
1 



is written as 



x = — J 2 /i|no sin 8 cos 8\L(f), 



where L(f) is the form factor given by 

I/-/- 1 ! 



L(f) 



l + P 



(67) 
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This expression shows that the number of spin reversals 
is proportional to the coefficient of friction /i, the axis in- 
clination <5, and the initial spin uq. Therefore, when the 
friction is large (close to the no-slip case) and the asym- 
metry between the shape axes and the principal inertia 
axes is large, the number of spin reversal becomes large. 

Here, I discuss the form factor L(f). When the long 
axis a is fixed, which responds to a fixed J2, enlarging 
/ correspond to reduction in J\ from Eq. (|44[) . Fur- 
thermore, from Eq. (|25|) . reducing J\ corresponds to 
a reduction in 6, which indicates that the rattleback be- 
comes long and slender. In contrast, reducing / results in 
the rattleback becoming more circular, such as J\ = Ji- 

The form factor L(f) has maximum at / = V^ + VE. 
Thus, it is observed that if the rattleback long and slen- 
der, rattle vibration about the short axis does not easily 
occur and the number of spin reversals is reduced. This 
also occurs, if the rattleback is more circular. The rela- 
tionship between a and b corresponding to the maximum 
of / is given by 



6 



75 



- 6. 



(68) 



For example, if a = 10, b becomes approximately 4.4. 
Fig. |3] shows the change in the form factor L(f) when 
a = 10 is fixed and b is varied from 1 to 10. It is observed 
that there is a maximum value around b ~ 4.4. 



A. Validity of the approximation i)o 3> v p 

In section Hill the approximation Vq 3> v p was used to 
derive effective equations of motion. 

The initial conditions are xq = yo = 0.01, xq = yo = 
and no = 0.05 with parameters a = 10, 6 = 3, and 
S = 0.03. Fig. [4] shows the behavior of vqi and v p i for 
jj, = 100, and Fig. O shows the behavior of V02 and v P 2 
for fi = 100. From Figs. |4]and[5l the comparison of the 




FIG. 4. Time evolution of vqi and i) p i for \i = 100 




FIG. 5. Time evolution of V02 and v P 2 for fj, = 100 



L(f(b)) 

0.35 




10 11 



FIG. 3. Form factor L(f) versus b for a = 10 



VI. NUMERICAL RESULTS 

@In this section, the results of computations based on 
the exact system described in Eqs. (|6|)- ([T0|) are presented 
and these results are compared with those based on the 
effective equations of motion described in Eqs. (|35 ]) -(|i? ]l . 

For the numerical simulations, the NDsolve command 
in M athematica is used. 



amplitudes shows that v v i ~ Vqi x 4q and v P 2 ~ V02 x m ; 
therefore, it can be safely assumed that vqi 3> v p i and 
V02 >• v P 2- Figs. [6] and [7] are the same as Figs. |4]and[5l 
respectively, but for \i = 200. In this case, it is similarly 




FIG. 6. Time evolution of i)oi and i) p i for /1 = 200 




FIG. 7. Time evolution of V02 and v P 2 for /i = 200 



shown that v p \ ~ Uoi x 155 anc ^ ^'p 2 ~ ^ 02 x TUB' therefore, 
vqi ^> v p i and Vq2 S> v P 2- Note that the approximation 
improves when [i increases. Thus, it can be concluded 
that the approximation Vq 3> v p is valid. 
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(b) 



FIG. 8. Time evolution of n(t) with (a) positive initial spin 
with n = 0.0015 (dotted line), 0.003 (dashed line), and 0.005 
(solid line) and (b) negative initial spin with no = —0.006 
(dotted line), -0.015 (dashed line) and -0.025 (solid line). 



B. n c i±, n C 2±, and spin reversal behavior 



In subsections IV Al and IV Bl of section [V] the critical 
values of spin necessary to increase rattle vibration and 
cause spin reversal are obtained. Here, computational 
results for these values and the spin reversal behavior 
which are obtained from a simulation based on the exact 
system are presented. 

The initial conditions are xq — yo = 0.01, and xq — 
yo = with parameters a = 10, b = 3, and 6 = 0.03. 
The coefficient of friction is set to fx = 75. For these ini- 
tial conditions, the following critical spin are obtained: 
n c2+ = 0.00358, n cl+ = 0.00179, n c i_ = -0.01, and 
n C 2- = —0.0201. Fig. [5] shows the computed evolution 
of n(t). In Fig. [5] (a), as discussed in subsection IV Al 
the dotted line with the initial spin hq = 0.0015 < n cl+ 
shows that rattle vibration does not increase. As dis- 
cussed in subsection IV Bl the dashed line with the initial 
spin n c i + < hq — 0.003 < n C 2+ shows that rattle vibra- 
tion increases but spin reversal does not occur. The solid 
line with the initial spin n C 2+ < uq = 0.005 shows that 
the spin reversal behavior agrees with that discussed in 
subsection IV Bl Fig. [8](b) shows the behavior of n(t) for 
negative initial spins uq = —0.006 (dotted line), —0.015 
(dashed line), and —0.025 (solid line). Even if the initial 
spin is negative, the same behavior occurs about rt c i_ 
and n C 2-- Thus, it is observed that spin reversal behav- 
iors are depended on these critical values. 



C. Number of spin reversals n r and friction 
coefficient /j, 

This subsection discusses how the number of spin re- 
versals n r changes depending on the friction coefficient fi. 
This number of spin reversals is obtained from a simula- 
tion based on the exact system described in Eqs. ©- ([TU| . 

The initial conditions are xq = yo = 0.01, ±q = yo = 0, 
and no = 0.05 with parameters a = 10, b — 3, and 
5 = 0.03. Fig. [9] (a) shows the behavior of n(t) for /i = 60 
(solid line), 130 (dotted line) and 200 (dashed line). It 
is observed that for each of these fj, values, the number 
of spin reversals increases to three, six, and nine times, 
respectively. Fig. [9] (b) shows the behavior of spin n(t) 
for fi = 500 (solid line), 1000 (dotted line), 1500 (dashed 
line). Note that the spin behavior is similar to that for 
the no-slip case as the friction increases. 



0.0.1 
0.02 



-0.02 
-0.04 
-0.06 



H=60 n=130 n=200 




ii y jj.lll'H; liliKi ; f Ami HmhV]...U( H lb 



H=500 n=1000 n=1500 




FIG. 9. Time evolution of n(t). (a) The initial values of n are 
60 (solid line), 130 (dotted line), and 200 (dashed line), (b) 
The initial values of /i are 500 (solid line), 1000 (dotted line), 
and 1500 (dashed line). 



We usually consider that when friction enlarge, energy 
loss become large, so that the reverse number decrease. 
The numerical simulations show the opposite thing to 
this intuition. In reality, the contact with the horizon- 
tal plane is not a point but an area, spinning torque by 
rotation n seems to effect the dynamics of rattleback as 
discussed by Garcia and HubbardQ. The adopted model 
does not include this spinning torque, so that energy loss 
depends only on the velocity v p which becomes small 
when friction coefficient \x becomes large, thus, it may 
be considered that the number of spin reversal increase. 
Moreover, in reality, when the spin reversal does not oc- 
cur, the contact point x p settles down to x = y = and 
z = 1, and the slip velocity v p is equal to zero. Then, 
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the rattleback stops spinning after a while. In contrast, in 
the numerical simulation, if we set the conditions v p = 0, 
x — x = x = 0, and y = y = y = 0, h becomes zero from 
Eq. (|26p. Therefore, the spin n(t) becomes constant and 
does not stop after the spin reversal ends, as shown in 
Fig. |9](a). It is considered that this phenomena are also 
depended on not considering spinning torque due to the 
frictional force by spinning. 

D. The approximate number of spin reversals n r:ap 
versus the exact number of spin reversals n r:ex 



the effective equations of motion is qualitatively correct. 
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In this subsection, the approximate number of spin re- 
versals n r:a p obtained from Eqs. (|65p and (|66[) in sub- 
section IV CI is compared with the exact number of spin 
reversals n r:ex obtained from a simulation based on the 
exact system described in Eqs. ([6])- (|10|) . 

The initial conditions are xo = yo = 0.01, xq = yo = 
and no — 0.05 with parameters a — 10, b — 3, and 5 = 
0.03. Fig. [10] shows the number of spin reversals n r for 
a range of values of fj, from 20 to 200. Crosses represent 
n r:ex and the boxes represent n r:ap . Up to fi ~ 70, both 




FIG. 10. The number of spin reversals n r versus \x: □ repre- 
sents n r -ap and x represents n r:ex 



values are almost identical, but, the difference increases 
as /i increases. It is assumed that the terms neglected in 
the approximation affect energy dissipation. Thus, for a 



given value of ji, 



l r:ap 



is larger than n r 



FIG. 11. Number of spin reversals n r - ex versus b 



VII. DISCUSSION AND CONCLUSIONS 

I have examined the behavior of a rattleback with vis- 
cous friction to answer the questions mentioned in the 
introduction. The following results were determined an- 
alytically and confirmed numerically. 

Critical value of the initial spin exist:n c i-|- and n C 2±- 
When the initial spin n is in the region ri cl _ < n < 
n c i +1 rattle vibration does not increase. When the ini- 
tial spin no is in the region n C 2- < tiq < n c i_ or 
n c i+ < riQ < n C 2+, rattle vibration increases but spin 
reversal does not occur. A numerical simulation based 
on the exact equations of motion confirmed the existence 
of these values. 

I theoretically obtained the number of spin reversals n r 
as a function of the coefficient of friction /i and the form 
factor L(f), which contains the ratio of a to b. From the 
expression of n r , it was found that the number of spin 
reversals increases as /x increases and a certain value of 
the ratio of a to & gives the maximum number of spin 
reversals. 

In this paper, viscous friction is adopted to perform 
a first examination of spin reversal. However, in reality 
it is assumed that spin reversal is associated with the 
effect of Coulomb friction. Furthermore, it seems that 
rolling friction and friction by spinning are also effective. 
Therefore, these frictions should be adopted to further 
understand the rattleback behavior. This is a subject for 
future investigations. 



E. Form factor L(f) and n r:ex 

This subsection discusses the relationship between the 
form factor defined in Eq. (|67[) and n r:ex . 

The initial conditions are xo — yo — 0.01, Xq — yo — 
0, and no = 0.05 with parameters a = 10, and 6 — 
0.03. Fig. [TTJ shows the behavior of n r:ex as a function 
of the parameter b for fj, = 100 and 200. Note that the b 
dependence of n r:ex is similar to that of the form factor 
presented in Fig. [3l Therefore, the estimate derived from 
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